Spherical Gravitational Collapse of Annihilating Dark Matter and the Minimum Mass 

of CDM Black Holes 



R. Ali Vanderveld and Ira Wasserman 
Center for Radiophysics and Space Research, Cornell University, Rhaca NY 14853 

Spherical gravitational collapse of a cold gas of annihilating particles involves a competition 
between two density dependent rates: the free-fall rate oc y/p and the (s-wave) annihilation rate 
oc p. Thus, there is a critical density p ann above which annihilation proceeds faster than free fall. 
Gravitational collapse of a cloud of (initial) mass M to a black hole is only possible if pBH = 
3M/4tt(2GM) 3 = 3/32ttG 3 M 2 < p ann , which implies that M > Af ann = (3/327rG 3 p an n) 1/2 . We 
compute the collapse and annihilation of cold, spherical gas clouds with either a homogeneous 
density profile or a specific radially-dependent profile that would lead to self-similar collapse if 
' there were no annihilation. These calculations verify that there exists a lower bound to the range 

of possible black hole masses that can form in cold dark matter (CDM) collapse. For a particle 
mass m and freeze-out temperature T; — m/xf, the minimum black hole mass is MaM « 1O 1O M X 
£Nj ■ (x/y^7/100ficDMft 2 3*s m (Gev)) , where g+s and <?* are degeneracy factors. Moreover, the formation 

of a black hole is accompanied by the annihilation of about M anrl of the original mass. Most of the 
annihilated mass is released in a burst lasting a time ~ GMbb, where Mbh is the black hole mass 
when collapse first occurs. A total annihilation luminosity ~ 10 59 erg s _1 could easily be reached, 
lasting a time ~ 100 seconds. The absence of astronomical observations of such spectacular events 
suggests either: (i) the branching ratio for CDM annihilation to e + e~ pairs or quarks < 10~ 10 , 
' while the branching ratio to vv is < 10 -5 ; or (ii) CDM is not made of annihilating particles, but 

may be in some non-annihilating form, such as axions; or (iii) CDM black holes never form. 
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| Non-baryonic matter dominates the matter content of the Universe today, with a density parameter f^cDM" ~ 0.15; 
the rest of the energy density is primarily in the form of dark energy In the favored cold dark matter (CDM) 

""^5 ' scenario based on weakly interacting massive particles (WIMPS), the dark matter is in the form of a gas of particles 
that is very cold today, with an abundance that was lowered substantially by annihilation when the temperature of 
I 1 the Universe was > m/xf, for a CDM particle mass m freezing out at temperature Tf — m/xf ~ 0.1m P, |j, [a] ■ The 
ACDM picture provides acceptable fits to the temperature and E-polarization fluctuations of the cosmic microwave 
background (CMB) radiation as well as to the large scale structure of the Universe On smaller scales, CDM 

halos seem to form cuspy cores jy, IToj ; dark matter annihilation has been proposed as one means of smoothing out 
• • . the cusps 

[HG3. CDM clumps are found to form in large scale structure simulations and to agglomerate to form 
halos ^3 ! annihilation within these clumps would lead to an enhanced flux of cosmic ray products at Earth 0, ^| . 

In this paper, we investigate a related question: How does annihilation alter spherically symmetric collapse of 
CDM? We consider, for the most part, the collapse of CDM clouds with zero or very small internal energy. In the 
absence of annihilation, the inevitable result would be complete gravitational collapse to a black hole. 

However, to form a black hole, a collapsing sphere of mass M must reach a mass density of order pbh = 
3M/47r(2GM) 3 = 3/32ttG 3 M 2 . When (s-wave) annihilation is included, the rate of CDM annihilation at a den- 
sity p is 

_ P(<TV) ann m 
^ ann — > I J- J 

m 

where (<7v) ann is the annihilation rate coefficient and m is the CDM particle mass. For comparison, the free fall rate 
at high densities is approximately 
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then r ann = H at a density 
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For p < Pannj free fall proceeds faster than annihilation, and vice- versa. Therefore, black hole formation is only 
possible for a mass M provided that pbh ^5 Paim or 

M>SFeM». (4) 

The annihilation of CDM particles in the early Universe leads to a residual density parameter ^cdm^ 2 provided that 

W—MxlO-O.V^( 1(Bn< Ji ' ), (5) 



where g*s and g* are degeneracy factors for the entropy and mass density when the temperature of the Universe is 
~ Tf, and 

x f re InC - iln(lnC) , (6) 

where C — 0.038(g/^/g^)m(av) ann /y/G, with g being the degeneracy factor for the CDM particles alone 0. Substi- 
tuting for (<7v) ann in equation @, we find 

^ 10 10 ^e f x f \ „ 1O 1O M / x f \ 

ann ~ m(GeV) U00Ocdm/i 2 9*s/V9;; ~ m(GeV) \Vz>g*s I yfc ) ' 



where m = m(GeV) GeV; furthermore, in our last result, we assume that ficDM^ 2 ~ 0.15, as has been determined 
from analyses of the CMB and of large scale structure. Our estimate of M ann depends on the CDM particle mass 
in several ways, since M ann oc Xf^fgljmg^s- In addition to the explicit inverse proportionality with m, there is a 
logarithmic dependence from Xf and a piecewise flat dependence via g+s and g+: for example, if m = 50 GeV, we 
would have enough particle species for g± ps g±s ~ 50, and M ann pa 4x 10 7 M Q (taking Xf ps 20 for this case). Plausibly, 
the value of the minimum black hole mass for CDM is comparable to the mass range deduced for black holes in the 
centers of galaxies 0] . It is remarkable (although perhaps fortuitous) that the value of M ann is even remotely close 
to an astronomically significant value. 

Spherical collapse is, of course, highly idealized. More realistically, we would expect a collapsing cloud to have 
substantial substructure as well as large scale shear flows induced by the tidal field of nearby fluctuations. However, 
we would still expect that CDM clouds that are, however improbably, nearly spherical and smooth would have the best 
chance to collapse to black holes. Thus, our lower bound on the CDM black hole mass ought to be robust. Moreover, 
shear enhances the collapse of a cloud (as is evident from the Raychauduri equation, e.g. llSl ') and anisotropic 
collapse of a cloud with internal fluctuations could trigger substantial fragmentation (e.g. []JJ,|23 and refs. therein). 
It is possible that the calculations we present apply to late time collapse within a larger cloud that has fragmented 
considerably. 

Section II presents calculations of the collapse of homogenous CDM spheres. Section III computes the collapse of 
a specific model with radial density dependence. In each case, we quantify when black holes may form, and when 
the collapsing mass will instead annihilate away. The detailed calculations are Newtonian, and assume instantaneous 
escape of the (relativistic) products of CDM annihilation. We expect that our main results for the critical mass for 
black hole formation, and for the mass-energy and luminosity released by annihilation, to remain valid within factors 
~ 1 in a general relativity calculation that also accounts for the opacity of the infalling CDM to the annihilation 
products. In Section IV, we discuss our results, including a brief account of what might happen in the collapse of 
CDM with predominantly p-wave annihilation. 

II. THE HOMOGENEOUS CASE 
A. Equations of motion 

Let us consider a uniform sphere of matter with time-dependent radius and mass, called R(t) and M(t), respectively; 
hence, the density at a time t is p{t) — 3M(t) /47ri? 3 (£). We assume that the homogeneous sphere collapses under its 
own self gravity, but we also include pressure forces. To maintain homogeneity, we assume that 



P(r,t)=P (t) 



1 R 2 (t) 



(8) 
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at time t and Eulerian radius r; we also assume that Po(t) = ifp 7 (i), where K is independent of time. The Euler 
equation then becomes 

- 2K(iM/^r^ GM 

R = - W - (9) 

If 7 > 4/3, then there exists a radius Rza at which R — 0. For a nonrclativistic gas we expect to have 7 = 5/3, so in 
this case we have the zero acceleration radius 

2K ( 3 x 2/3 



RzA= GM^{^) ' (10) 

This is a minimum radius at which the pressure force becomes as great as the gravitational force such that the collapse 
will "bounce" back outward. In the absence of annihilation, we also have the conserved energy 

1 - 9 GM /3M\ 2/3 K , . 
E=-R 2 + — . 11 

2 R \ 4tt J R? V 1 

From this we see that the cloud will bounce at a radius R = \Rzk in the absence of annihilation. 

We note that there may be three possible physical sources of "pressure": (i) internal motions associated with 
small scale density inhomogeneities within a cloud; (ii) residual but small entropy of the CDM particles; (iii) trapped 
annihilation products. 

Of these three sources, the first is potentially the most important, at least on large scales, but, as discussed in the 
introduction, we are focusing on fragments that have the best chance for collapse. Thus, we assume that our collapsing 
cloud is smooth; the larger scale motions associated with inhomogeneities of the parent system are assumed to be in 
the form of orbital motion of the subclumps into which it has fragmented by the time our smooth cloud collapses to 
very high density. 

The residual entropy of the CDM is small and depends on the scattering cross section of thermalized particles in the 
early Universe with the CDM. Assuming a cross-section ft(o'v) Sinn (T/m) 2 , the rate of energy transfer from thermal 
particles to CDM is U t = C,{'i)gtft{o'v) a _ nn T 7 /T^rn 3 per CDM particle, where gt is the spin weight of relativistic 
thermal particles that scatter from the CDM and ft is a numerical factor of order one. The CDM therefore remains 
in "kinetic" thermal contact until T t = 2U t /3T = H = (8tt 3 ^GT 4 /90) 1/2 , which implies kinetic freeze-out at 

Lgg^gVV^ _ 3 g^KGeVjT f lOOiW V" . .... 

K ~WW1^F' WW {—^—) • (12) 

the value of K associated with this residual thermal energy is then 



m[p(T x )]2/3 



4.9 x 10 8 GeV- 8/3 9t ft „,„ ( 



!/4 f l/4 / ^. \ V4 / q 15 \ 2/3 



[m(GeV)] 7 /4 5 VV^ VWO^cdm/i 2 / V^cdm^ 2 



(13) 



where g+sn is the degeneracy factor for the entropy at temperature Tk- To put this in context, the CDM phase space 
density is - if^m^ 4 - Kr 13 [77i(GeV)]- n / 8 . The CDM phase space density and, therefore, -Kcdm are conserved 
as long as the CDM particles are collisionless. 

The third source of pressure would be from trapped products of the annihilating dark matter. For the most part, 
we shall assume that these annihilation products escape without being trapped, and this further assumption will be 
discussed below. (We shall also re-examine this assumption semi-quantitatively in our brief discussion of the collapse 
of CDM that annihilates predominantly via p-wave interactions.) Any products that are trapped should be extremely 
relativistic and therefore would form a 7 = 4/3 gas. At fixed K, a 7 = 4/3 gas either stops the collapse at all times, 
or not at all. 

Annihilation implies a mass loss rate 

M = - 3 <7>-f. (14) 

Here, we assume that the annihilation products escape the collapsing cloud completely. This and equation 10 give us 
our system of coupled differential equations which we must solve. These look much simpler once nondimensionalized 
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with a change to new mass, radius and time units. Let the mass unit be the initial mass Mq such that the new mass 
variable is M(t) = M(t)/Al . Then introduce a radius unit R u and a time unit y/fi^/GMo for the new variables R 
and t, respectively. Then the above equations become 



d 2 R _ R ZA M 2 / 3 M 
~dr 2 ~ R U R 3 # 



(15) 



dM 3(av) ann IGMoM 2 ( R ann \ 3/2 M 2 

(16) 



dr AirmG y R^ E? \ Ru J R? ' 

here the characteristic radius at which annihilation becomes important is defined to be 



Rann — 



3(0-w)anii f Mq X ' 



^/ 3 irGm \M ann 

3.8 x 10 15 cm/ x f \ ( M Q \ 1/3 



m(GeV) \lOQn CBM h 2 g+ s /^g;J \M 



(17) 



We choose R u — \J R 2 nn + R\^ letting k — Rza/Ru and i? an n = R U V^ — k 2 , the equations describing the collapse 
become 

d 2 R kM 2 ' 3 M 
'dr 2 ~ ~R3 B 2 



dM (1 - k 2 ) 3 ^M 



2 



dr R3 
The characteristic timescale for the cloud's evolution is 



(18) 



t u 



R 3 3(crw) an „ f R u 



GMq 4nGm \R a 
2.0 x 10 5 sec 



3/2 



x f 



m(Gev)(l- K 2 )3/ 4 \l00n C BMh 2 g, S /V^ 
For cases where pressure forces are unimportant, n — » 0; for K = Kcbm 

i?zA _ k „ 1.3 x lQ-^V^gff 12 /10WW^\ 13/12 / 0T5 \ 2/3 (M m N 2/3 



i? a „n v/1 -k 2 g% 3 K gl /3 [m(GeV)} 5 / 12 \ x f J V^cdm/i 2 / V M < 



(19) 



(20) 



which is indeed very small. 

The solution of equations (|18fl depends on the value of n and the initial conditions. Well before annihilation becomes 
important, M remains constant, and 

f^y=— +i+e. (2D 

\dr J B 2 R 

For E = 0, we can integrate to find 

2 /- K\ 3 /2 /. K Nl/2 



V2\r-T bounce \ = -[R--J +K { R -2J ' (22) 



where the bounce happens when t = Tb ounce . We use equation (1221) to estimate the initial time for our numerical 
solutions. 

As for the initial conditions, M(t) = 1 at the start and R(t) should ideally fall from a very large distance, where 
dR/dr w 0. However, since annihilation occurs only at small radii R < 1, we can initialize numerical calculations 

to excellent accuracy by setting Mq — 1 as long as the initial value Rq is sufficiently large. We can then find the 
corresponding initial infall velocity by assuming previously negligible annihilation; in this limit, we can use the above 
energy conservation equation to find 




(23) 



which will be used in the numerical work described below. 
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B. Numerical calculation 



Equations (|18fl were solved using a fourth-order Runge-Kutta routine |2l|. The initial conditions for the mass, 
radius and velocity were set as discussed above; Ro = 10 turned out to be sufficiently large. The results, for different 
k values, are displayed in Figure Q Note how the minimum R and the final M values both appear to scale with 
k. Also, we see that the k = case yields a power-law relationship between M and R for small M and R. We can 
understand these relationships analytically by looking at the small k and small M limit of the equations of motion. 




1.179 1.18 1.181 1.182 

log(r) 




log(R) 



FIG. 1: The top two plots show log(i?) and log(M) versus log(r), for a time interval near the bounce and for k = 0.03, 0.01, 0.003 
and 0. The bottom plot is log(M) versus log(/j) for the same four values of k. The initial radius was Ro = 10. 

Let us change the independent variable in our equations from r to x = X/M while defining the new dependent 
variable y = R~ 2 . Rearranging the equations of motion yields 

This form is more useful for finding the behavior of the solution when M — » 0, which is equivalent to x — > oo. In this 
case, let us look for a solution of the form 

y = Ax + Bx a + ... (25) 

where the solution will only be sensible if the second term is much smaller than the first. By inserting this into 
equation we find that 




y ~Ax-8J- + ... (26) 



for large x; the coefficient A is not determined by this argument. Note that, by returning to the original variables, 
this implies that M/R 2 approaches a constant as the sphere collapses. Putting k back into the equations, but still 
requiring it to be small, we can take the R — > limit in order to predict a final mass Mf cx k 3 and a minimum radius 
\Rza <x k 3 / 2 . These scalings agree with the numerical solutions. 

The asymptotic scaling M j R 2 — > constant for an annihilating cloud implies that dR/dt — > constant at late times. 
Before annihilation becomes dominant, the cloud collapses in free fall; this phase ends when the cloud attains a density 
~ Pann- At that time, the cloud collapses at a rate dR/dt w u ann = (Mo/Mann) 1 / 3 ; subsequently, the cloud loses 
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mass, its self-acceleration diminishes, and it continues to collapse at roughly this speed. In doing so, its density is 
determined by the condition that the rates of annihilation and collapse remain about the same, i.e. 



PWann 3M(av) 



m AirFPm R 

which implies 

M 47rm« ann 



R 2 3(av) 



(27) 



(28) 



Equation corresponds to M/R 2 = A — 4 1/3 . 

At this point, we pause to consider the circumstances under which our assumption that annihilation products escape 
promptly is justified. Let us suppose that, when first produced with individual particle energy m, the relativistic 
annihilation products interact with the CDM particles with a cross section f s (o'v} ann , where we expect / s ~ 1. Then 
the optical depth of the collapsing cloud to the relativistic annihilation products is r = f s p{av) ann R/m. The density 
at which the cloud first becomes opaque to these products is 

f _3/ 2 /^annV /2 ,-3/2 ( M V'* (9Q , 

Pl=Pa„n/ s (^-) =PBH/ S l^-J - (29) 

Thus, we see that the cloud remains transparent to the annihilation products for cloud masses well above and well 
below M ann . For cloud masses » M ann , the optical depth rises to t ~ 1 just as the cloud reaches p ann ~ Pbh- For 
these cases, the collapse of the cloud could be altered slightly, delaying the final outcome - total collapse or total 
annihilation - for a short time, but probably not for longer than the characteristic free fall or annihilation time. 



C. Black hole formation 



In order to analyze the black hole formation process, define the dimensionless "horizon parameter" 

Kr) - ^ - , (30) 

1 1 R( r ) Rij) M ' v ' 

which increases during the infall phase, reaches a maximum at the bounce and then decreases during the outflow 
phase. For k^O, this simplifies to 

h( x 2^/3 \ 2/3 M(t) 

Kt) = ^\-^mI) ~R{t) ' (31) 

and for K = 0we get 

V ' WmVGMoJ R(t) v ' 

We can easily find the maximum value of h(r) for any particular value of n, which ranges from h max {K = 0) ~ 0.8 to 
h m ax(n = 1) = 2. This is shown in Figure|21 In order for a black hole to form, we would require that 

T > 1 < 33 > 

and this corresponds to a restriction on Mo. For n ^ 0, we require 

T/ K \ 3 ' 4 



for black hole formation, and similarly for k = we would require 

3 M ann / 1 \ 3/2 y/2M, 



ti t ^ " \" "/ aim / * I _ v ~^~jmn Ci^) 

°~ inmG 2 \2h max (0)J ~ [h max (0)f 2 ' 1 j 

Therefore there is a definite minimum initial mass ~ (1 — 2)M ann required in order for the self-gravity of the system 
to overcome the annihilation and pressure. 
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FIG. 2: The top plot shows h ma x versus k and the bottom plot shows the values of M and R at this value of h, also versus k. 
Initial radius Ro = 10. 



III. A RADIAL DENSITY PROFILE 



A. Equations of motion 



A purely collisionless gas cloud, defined to have zero pressure and to be non-annihilating, collapsing under its own 
self-gravity can be visualized as a set of mass shells in free fall. We can label each shell by the mass it encloses, m s , 
such that its radius at time t is r(m s ,t). As long as no two shells cross, m s remains unchanged and can be used 
as a comoving coordinate; this remains true when annihilation is included if we choose m s to be the mass enclosed 
long before annihilation becomes important. For purely collisionless infall, the motion of an individual mass shell is 
determined by its conserved energy E(m s ) and its "bang time" to(m s ), which is defined by r(m.^tp(m s )) = 0. 

Self-similar solutions for the collapse of a collisionless cloud have been constructed by Penston [22] and Lynden-Bcll 
and Lemos j^, hereafter LL. These models correspond to E{m s ) = and a power-law to(m s ) = <o + Bm b ; the exact 
solution is then 



r(m s ,t) = 



/9Gr 



V 2 



1/3 



(t (m s ) - 1) 



_ ^2/3 = f 9Gm s 



1/3 



{Bm b s +t -ty 



(36) 



The density and velocity profiles corresponding to this solution are 

1 



p(m s ,t) 



v(m s ,t) 



1 



47r[r(m s , t)] 2 (dr{m s ,t)/dm) 6irG(Bm b +t - t)[Bm b (l + 2b) +t -t] 

1/3 



2Gr 



r(m s ,t) 



4Gm, 



3(Bm b +t -t) 



(37) 



where the density profile differs from LL equation (2.9), which has a typo. Equations (|36|l and (|37|l describe a collapse 
profile consisting of a shrinking core of mass 



and radius 



m c ~ [(t -t)/B]V h 



r c ~ (18G) 1 /3 B -V36( to _ t )(l+26)/36 



(38) 



(39) 



Within the core, the density is nearly uniform, p ~ p c = l/67rG(£o — t) 2 , and the velocity profile is nearly linear, 
v(r, t) ~ —r^/WKGp~ c . This is surrounded by a halo, in which the density and velocity profiles are approximately 
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powerlaws: 

p ~ p c (m c /m s ) 2b ~ p c (r c /r) 6b /( 1+2b ) 

» ~ DcWm)^ 3 ~ WcCrc/r)^- 1 )/^ 4 - 26 ) . (40) 

In the LL solution, the core density rises to infinity as t — > to, so black hole formation is inevitable in this model. 

Our goal is to numerically compute the solution to this problem in the presence of annihilation. As in the preceding 
section, we will choose initial conditions such that annihilation will have been negligible previously, making our use 
of the LL solution a physically meaningful choice. The above solution has its largest density at the center, so we 
must specify the central density to be small enough such that it is sufficiently below the critical density specified in 
equation The central density at the initial time t = is 



f 



For our numerical simulations, we require that 



Pann = -^P(m s = 0,t = 0) J (42) 



this puts a constraint on the constant t$: 



6tt Gm 



(43) 



where we choose the constant factor a such that the initial central density is sufficiently small. We chose a = 3/47rx I0 6 , 
which starts the collapse at such an early time that the central density is approximately a million times smaller than 
Panm and therefore the collapse is still in freefall with utterly negligible annihilation. In this way, we will be justified 
in initializing our simulation with data from the LL solution. 

In order to determine the other constant B for our initial conditions, we can now specify an initial density contrast 
between the center and the outer edge of the gas cloud; specifying this contrast is also necessary for numerical purposes 
because the LL solution extends off to infinity. Let the central density be a factor of fa larger than that of the edge 
at the initial time; this assumption, along with density equation 1)3 7|). tells us that we must have 




(44) 

where we have once again chosen to call the total initial mass Mq and 7 is the constant 



y/(b + l) 2 + (2b+l)(ft - 1)- (6+1) . . 

7 26+1 ' 1 ' 

The above conditions give us the constants and initial data necessary for a numerical computation. 

Although we are using the initial contrast fa for our initial data, the physically important density contrast is what 
we will call (3 a nm which is the contrast when the center begins annihilating. This is the parameter which will become 
important when discussing black hole formation for this system. We can write fa as a function of P ann and a. First, 
define the annihilation time t ann to be the time at which the central density would reach the annihilation density in 
the LL model: 

1 

67rG(£ - tannY 

Then (3 ann = p ann /p(M , t ann ), which is 



fa nn = [BMq y/6nG Pann + 1 (1 + 2b) BM$ y/6nG Pann + 1 (47) 



and we can solve this equation for BM^y/6nGp a nn to get 




S « \/i^ = -T^ + W^+ ' < 48 > 
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The initial density contrast is 



Pi 



( ^ + i)(li± ^ + i) 



which can be rewritten as 

ft = 1 + 2 (1 + 6) 



(BM^y/^Up a 



V toV^Gp ann 

plugging in what we know from the equations above, we find 



+ (1 + 24 +™M 



anti 
ann 



(49) 



(50) 



ft = l + 2 ( l + 6 h / — 



' ftann 1 

1 + 26 



1 + 6 
1 + 26 



1 



1 + 26 



+ (1 + 26) 



—) 

Sua J 



'Pa 



1 



1 + 26 



+ 



1 + 26 



1 



1 + 26 



(51) 



We will use ft and a in the calculation of initial conditions, but whether or not a black hole forms will be sensitive 
only to our choice of [3 ann . 

The addition of annihilation requires that we add the mass loss equation 



dp(m s ,t) _ (av)z 



dt 



m 



-p(m a ,t)fj,(m s ,t) 



(52) 



where the comoving variable m s now must denote the mass enclosed at the start of the collapse, while the new 
dimensionless variable p(m s ,t) is the fraction of mass at a coordinate m s which has not yet annihilated at the later 
time t. Then the density p is simply 



p{m s ,t) 



p(m s ,t) 



Airr 2 (m s , t)dr(m s , t)/dm 
In the presence of annihilation, the collapse equation becomes 

d 2 r(m s ,t) _ G 



dt 2 



r 2 (m s ,t) 



/ dm' s 
Jo 



H{m' 8 ,t) . 



(53) 



(54) 



B. Numerical calculation 



To solve our system of coupled partial differential equations numerically, we divided the collapsing sphere into N 
shells, initially with equal masses and with the nth such shell having a time-dependent mass Am n (t) and radius r n (t). 
If n > 1 , then the spatially discretized equations of motion for this shell are 



d 2 r n {t) 
dt 2 



G 



r n (ty 



i=l 



dAm n (t) 3(cru) ann ( dm 

-Am n (t) 



dt 



4irm 

3(o-«)ann 
47T77J 



dr 3 
Am n (t) 2 



r„(i) 3 -r„_i(t) 3 



(55) 



where the central "shell" (n = 1) is treated as a homogenous sphere as in the previous section, and thus will have the 
following equations of motion: 



G 



:Ami(t) 



d 2 n{t) 

dt 2 ' ri (t) 2 

dAmi(t) 3(<rv) ann / Ami(t) 2 

dt 4irm \ n(t) 3 



(56) 
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Our requirement in the previous subsection that a be large also requires that this innermost shell is significantly 
larger than its annihilation radius at t = 0. For a — 3/4ir x 10 6 , ri(t = 0) is one hundred times greater than the 
annihilation radius for a homogenous sphere of mass Ami(i = 0). 

In order to nondimcnsionalize this system, we choose convenient units as before with the homogenous case; let the 
mass unit be the initial shell mass Amo = Mq/N and let the length unit be the previously-defined annihilation length 
for the central shell 



Ra 



We also choose a corresponding time unit 



3(av) 



AirmG 



2/3 



( GMp 
V N 



1/3 



R 3 



GAm 



(57) 



(58) 



such that the equations take on a simpler form in terms of the new dimensionless variables Am„ = Am„/Amo, 
f n = r n /R ann and r = t/t c : 



d 2 f n (r) 
dr 2 

dAm„(r) 
dr 



1 " 



Am n (r) 



r„(r) 3 - r n -i{T)'- : 



(59) 



for n > 1 and 



d 2 fi(r) _ Ami(r) 
dr 2 ~ 
G?Ami(r) 



dr 



ri(r) 2 
Ami(r) 2 
fi(r) 3 



(60) 



for n = 1. These equations were integrated in time using a fourth order Runge-Kutta routine similar to that in the 
above section, with the initial time t = 0, initial shell masses Am„ = 1 and initial radii given by equation (|36f) to be 



f n (r = 0) 



/ 47rcm^ 



b \ 2 / 3 
+ 1 



The initial infall velocities were chosen to be 



df„ 
dr 



A 



5> 



(61) 



(62) 



which corresponds to zero energy infall. 

As we would expect from the homogeneous case, each shell reached both zero size and zero mass at the same time 
in the absence of black hole formation, i.e. for clouds which annihilate completely. Some results for this simulation 
are shown in Figure for b = 2, a = x 10 6 , f3i = 100 and N = 1000. We chose to use one thousand shells 

for the simulation because this N provided sufficient spatial resolution while still allowing the code to run within a 
reasonable amount of time. 

It is important to note that, for this set of parameters, there appears to be a transition around a time r ~ 500. 
This is simply the annihilation time t ann , defined in equation 146|) to be the time at which the center reaches the 
annihilation density: 



1 



y/6-KGpa 




(63) 



using to = y /, a/6ir((av) ann /Gm), a = 3/Attx 10 6 and our nondimensional time variable, we find that T ann = t ann /t c w 
471. This is the same value as we see in the numerical solution, and therefore we can conclude that neglecting 
annihilation effects until reaching the critical length scale is a good approximation for this choice of parameters. 
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FIG. 3: Data for a simulation run in which a black hole does not form, with b — 2, N — 1000, a — 3/4-7T x 10 6 and (3 = 100. 
The top plot shows the log of the dimensionless density p versus the log of the radius f for several times r = 0, 200, 500, 1000 
and 2000. The bottom plot shows the dimensionless mass enclosed versus radius for the same r values. 



C. Black hole formation 



Our main goal is to generalize our investigation of black hole formation in uniform density collapse to non-uniform 
collapse. To that end, at each time step, we check to see if any shell satisfies the condition 



2Gr 



> 1 ; 
r 

in terms of our nondimensional variables and spatial discretization, this condition becomes 

1/3 



27V 2 M, 



TO,; > 1 



(64) 



(65) 



If this condition is true for a given shell, then we stop integrating the equations of motion for all shells inside it. In 
other words, we find the location of the event horizon and shut off the annihilation and collapse of the mass inside 
it. Then, if a black hole does form, there will be an initial seed hole that subsequently will accrete mass from the 
remaining (annihilating) matter outside. If a black hole does not form, then all of the shells will eventually annihilate 
away until no mass is left. Again, we shall sec that the initial cloud mass Mo must be above some critical minimum 
value, M c , in order for the cloud to collapse to a black hole. 

In general, the mass scale M c depends on f3 ann and b. To determine M c for a given choice of b and (5 a nm we increased 
our inputted value of the mass ratio Mo/M ann until a black hole could form. (Of course, (3 ann — 1 reproduced the 
results found for the homogenous case, independent of the value of b.) As we increased the contrast (3 ann , we found 
that this minimum mass, or M c , also increased via a rough power law. These results, for several b values, are plotted 
in Figure Note that the exponent of the power law depends on b. 

We can estimate the critical masses in Figure 0] by using the LL solution and simultaneously requiring that a black 
hole forms before the annihilation time t ann . In the LL solution, 



2Gm s 



AGm s 



_3(Bm b s +t -t) 



2/3 



(66) 
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FIG. 4: Simulation data plot of \og(M c /M ann ) versus log(/3 ann ) for b = 0.1, 1 and 2. 
and this has a maximum value when 



//;.,. 



tp-t 
5(6-1) 



1/6 



(67) 



note that this becomes singular as b — > 1, and therefore we will assume b > 1 for the time being. Inserting equation 
(|67|l into our equation for <f> gives us the maximum value 



4G (6-1) 



l-l/b 



2/3 



1 at the time of initial black hole formation, tsH, i-e. when 

b 

6-1 



t = t 



BH = <-0 



4GN ^ r 

T 



1 b 

Qb-l\}b-l 



(68) 



(69) 



In order to produce a black hole before annihilation becomes important, we want 4> m ax = 1 and t < t ann ; this happens 
if 



B < 



4G\ b (6-1)"- 1 



3 J b b (t - t ann ) 
Using equation (|48(l . we can rewrite this condition as 



6-1 



(70) 



M > 



36 (t Q - t a 



4G (6-1) 



1-1/6 



6+1 
26+1 



1/26 



(3a, 



1 



26+1 



6+1 
26+1 



1/6 



Mr. 



(71) 



For large f3 ann , this yields the power law M c oc /3, 

If 6 < 1, the situation is somewhat different. Note that, for any b value, 



d /2G 



dm s 



oc t - t + (1 - 6) Bm^ , 



(72) 
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which is positive for b < 1 and to — t > 0, indicating that <f> is an increasing function of m s under these circumstances. 
In that case, (j> is largest at the edge, where m s = M . Following the same line of reasoning as for the 6 > 1 case, 



requiring that 



1 and t < t a 



we find that 



3 (to t an n) 

4G 




0a 



1 



26+1 2b + 1 



(73) 



for the b < 1 case. Using the data in Figure we have found that this estimate is too large by a factor ~ 2 for the 

1/2 

b = 0.1 case, even though the data do follow the predicted power law M c oc f3 ann for larger values of j3 ann \ in general, 
the estimate of equation (|73|l is too large by factors of order unity for b values less than one. Hence, for b < 1, we can 
conclude that it is possible for black holes to form at a later time i = t ann + At, such that M c is reduced by a factor 

Of (to — tann ~ At)/ (to — t an n)- 

The simulation also tracked the masses of both the initial seed black holes and the end-state black holes that they 
eventually evolved to. Theoretically, for 6 < 1, we expect that the initial black hole will have a mass of order Ma, 
and this prediction was verified by the simulation. For b > 1 , we should find the initial black hole mass from equation 



M BH = 



to — tBH 

B(b-1) 



1/6 



(74) 



plugging in what we know for tsH gives us 



M BH = 



bM n 



' 6+1 
26+1 



Pann ~ 1 b + 1 



26+ 1 26 +1 



oc A L 



( Mo 



\M a 



(75) 



which agrees with the simulation data. 

Finding the final black hole mass is equivalent to finding the total mass change due to annihilation, which we will 
call AM. Suppose that a black hole forms with a mass that is Mbh initially. The mass density at this time is pbh ~ 
3/327rG 3 M^ H , so annihilation proceeds at a rate w pbh (c'u)ann Mb#/to ~ 3(av) ann /32TrG 3 mMBH ~ M ann /GMBH- 
Up to numerical factors, we therefore expect a peak annihilation luminosity — M ~ M ann /GMBH an d a total mass 
loss to annihilation AM ~ —M ann - Thus, whenever a black hole forms, a similar amount of energy is released in the 
form of annihilation products. 

To study this more completely, we consider collapsing clouds well above the critical mass for black hole formation. 
The easiest case is the collapse of a homogenous cloud (6 = 0). If we assume that the initial mass is significantly 
above M c , then we can assume that annihilation will not affect the trajectories of the shells as they fall inward. For 
the uniform case 



Pit) 



1 



67rG(t -t) 



(76) 



and a black hole forms when to — tBH = 4GMq/3. The total mass annihilated is 



AM = 



(crv) ann Mo 

GirGm 

(<™)ann _ 

'8irG 2 m ~ 



dt 



oo (to — ty 



(o^)annMo 

67rGm 



1 



*0 — tBH 



(77) 



Therefore, the total mass that is lost to annihilation in the homogenous case is of order M ann , for any initial mass 
that is greater than the minimum mass M c . 

All cases with 6 < 1 are similar to the 6 = case in that the event horizon first forms at the outer edge. In this 
case, annihilation ends when 



to — tBH 



4M 



BMq ; 



(78) 



using the LL solution for p(m s , t) we find 



AM 



WnbGmBM^ 1 Jo x b 



dx 



In 



(26 + 1) BMoX b - BMg + 4M /3 



BMgx b - BM b + 4M /3 



(79) 
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In the limit that 6 goes to zero, equation (|79fl yields the same answer as before: AM — > —2M ann /3. 

For 6 > 1, the situation is more complicated because annihilation continues after the initial black hole formation. 
For this case, it is useful to think in terms of the luminosity as a function of time, which is 



Tl .\ - dAM Wann , , t , 

L(t) = — — = / p(m s ,t)dm 

at m 



cry 



Mo 

M min {t) 

° dm* 

(80) 



67rGm J Mmm (t) ( Bm s + *o - *) [(26 + 1) Bm b + 1 - t] ' 

where M m i n (t) is the mass within the event horizon, which doesn't contribute to the luminosity. Let tsH be the time 
when a black hole first forms; then for t < tsn, we have M m i n (t) = 0, and for t > tsH we have 

AGM min (t) , ,,«, , s 
= B [M mm (t)] b +t -t. (81) 

For very early times, when t < tBH and to — t >> BMq, there is an initial phase where we can approximate 
/ dAM\ (av) ann f Uo dm 



( dAM \ Hang f M0 dm, 

V dt J early 6irGm J (Bm b s + 1 - t) [(26 + 1) Bm b s + 1 - t] 



M 



(to-t) 2 



6irm 

For smaller to — t, but still before black hole formation, we have a another phase where 
/ dAM\ (<jv) ann f M ° dm s 



(82) 



(83) 



V dt J t<tBH 6vGm{to-tf Jo [Bm b /(t -t) + l][(2b+l)Bm b /(t -t) + l]' 

if we define x b = Bm b s /(to — t) and assume that {Mo/MBH) b >> 1, then this becomes 

/ dAM \ b 2 M ann {to-t BH \ 2 - 1/b r x dx 

{~^r) t<tBH ^~2(b-lf- 1/b GM BH { h-t ) Jo (^ + 1)[(26+1)^ + 1] ' 

The dAM/dt of this second phase increases with time until it ends right before t = tg^, and so we get a peak 
luminosity of 

• ^ b 2 M ann [°° dx 

M peak „ 2 ^ _ i)2 _ 1/ft ^ + i} m + 1)xb + 1 y m 



(86) 



Right at t = tsH, we know that M m i n = Mbh, and so then we have 

( dAM \ b 2 M ann f 00 dx 

V dt J t=tBH * 2(b-l) 2 - 1/b GM BH J( b -i)-^ (* b + l) [(26 +1)^ + 1] ' 

and thus, in our simple model, dAM/dt drops by a factor of order one right after a black hole forms. 
The total mass that is annihilated before tsH, which we call AM\, is simply 

f tBH , dAM(t) 

AM 1 = / dt -li ; (87) 

J-oc dt 

the mass annihilated from the phase where dAM/dt w (dAM/dt) ear i y is approximately 



AM, 



early 



2bM ann [M BH ^ b 1 



Mo 

which is very small for Mq/Mbh >> 1, and therefore we will neglect this term; then equation l|87|) becomes 

f tBH , fdAM(t)\ 2b 2 M ann f 00 dx 
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where we took the upper mass limit to infinity, which is justified for Mq » Mbh and 6 > 1. 

After tsH, we need to use the more general equations (|8Uf) and JSTJ, but for very late times, when t — to >> 
GM m i n (t), we can approximate 

(t-t \ 1/b 

M mm {t)~ (— g 1 ) + . (9°) 



where 

= AM min {t) 
€ 3b(t-t a ) ' 

then, in this very late phase, we can define the integration variable x = m s /M m i n (t) to get 

( dAM \ {(Jv) & nnM mm {t) f 00 dx 

V dt J late ~ 6TrGm(t-t ) 2 h [(1 + e) x b - 1] [(26 + 1) x b - 1] 



AM 2 w 1 fd*,. I . 



(91) 



2(6- 1) 1_1/& GA/ BH V *-*o / \t -t B H 
The total mass annihilated after which we call AM2, is given by 

AM 2 = / dms / — — — — — , , (93) 

6nGm J Mbh J tBH (Bm b s + 1 - 1) [(2b + 1) Bm b s + 1 - 1] K ' 

where tf(m s ) is the time at which a mass shell with a coordinate m s reaches the event horizon: 

^pL = Bm b s +t -t f (m s )-, (94) 
if we use the integration variable x = m s /MBH, we find that 



(95) 



for M » M B h and 6 > 1. 

We can integrate the above equations numerically to get AMi, AM2 and M pea k, and we show the results for 
several values of 6 in Table [I] We also provide a sample lightcurve in Figure [5] which shows the scaled luminosity 



TABLE I: Data for several values of b 



b 


AMl/Mann 


AM 2 /M ann 


GMpeakMsH /M ann 


2 


-1.29 


-1.08 


-0.97 


3 


-1.01 


-0.57 


-0.76 


4 


-0.91 


-0.39 


-0.68 



LMbh /M ann as a function of the scaled time t(s)(10 7 M Q /MBH) ~ (50 s)(t/GMBH)- Our calculated light curve 
ignores the effects of opacity and general relativity. The former should be relatively unimportant when a black hole 
forms, as was discussed for the uniform case in Section IIB. (When 6 < 1, the uniform case is a good approximation, 
and when 6 > 1, most of the optical depth for annihilation products is near the hole.) General relativity would tend 
to lower the luminosity peak and smooth it out. Within factors of ~ 2 or 3, we still expect a lightcurve resembling 
Figure |31 and a total annihilated mass ~ M a 



L ann- 



IV. DISCUSSION 
A. P-Wave Annihilation 

P-wave annihilation has {av) ann — uo{v 2 ). In the early Universe, (v 2 ) — T/m, and following we find 

r-2 / X f 



a « 1.7 x 1(T 7 GeV- 2 — ± —= , (96) 
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FIG. 5: Plot of L(t) versus t for b = 2. 



where Xf w C — § ln(lnC) with C w 0.076(g/ v ^)moo/\/G- Assuming that thermal particles scatter from the CDM 
particles with a cross-section f t ao(T/m) 2 , we find that the CDM remains in kinetic equilibrium until 



T K « 1.6 x 10-3 GeV ^ [ ^ )] 

consequently 



y c 4 fm(GeV)l 3 / 4 / 10 3 Ocdm^ 2 



1/4 



5 ! /4 / t 1/4 V 4 



(97) 



Kcdm ~ 10 9 GcV~ 8 / 3 gt 1/4 2/3 ^ , 2 (0^5) ■ 08) 



As long as the collapse proceeds adiabatically, (w 2 ) = KcdmP 2 ^ 3 , and for pure p-wave annihilation we find 

^ 2.4 x 1Q»M / gt 3/28 /f /2 V 3/14 \ / 4 \ 15/28 f 0-15 ^ 2/7 
ann ~ [m(GeV)] 33 / 2 s ^ glT*9lsK J V 103Qcdm/i V V^cdm^ 2 . 



(99) 



For example, if m = 50 GeV and g t w 5* 5^5 w g*s/f w 50, we find M ann w 10 6 / t 3 ^ 28 M©. Not suprisingly, the 
value of M ann is considerably smaller than for s-wave annihilation, under the assumption of adiabatic collapse: p-wave 
annihilation is less efficient, and therefore the cloud must reach higher density to annihilate, which lowers the critical 
mass for black hole formation. 

In actuality, we do not expect the cloud to collapse adiabatically, and the critical mass is closer to the value for 
s-wave annihilation at (av) ann = <7o- The reason is that p ann is large enough that clouds become opaque to their 
annihilation products before they either annihilate or collapse catastrophically. Defining p\^d to be the density at 
which the optical depth is one in adiabatic collapse, we find 

P1M = Pan„/ s - 3/2 (^p 2 / 3 n) 3/2 (^) ^ = /Wf 8/8 M£ ) 3/2 (j^) ^ ! ( 100 ) 



since 



K n ys . 1-2 xlO- 3 (g?V^r\ (l^cnufA ^ ( 0-15 \ 2/7 

iWp ann ~ [m(GeV)]5/28 ( ;w II ^ ] ' ( 101 ) 
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we see that the clouds become opaque before reaching their final destinies. 

Once the clouds become opaque, they start to heat up. To understand what happens qualitatively, we introduce 
new density and mass scales 



Pi = 



Mi 



3/X 




32ttG 3 pi 16ttG 2 to 



(102) 



in terms of which we find 



P y ■ (M_^< ' 
Mi 



4 



Pi 

2GM 



R \pi 

we can eliminate the density from these equations to find 



Mi 



M 
Mi 



1/2 



(103) 



(104) 



Thus, we see that clouds with M > Mi collapse to black holes before becoming opaque, and vice-versa. We also see 
that during adiabatic collapse, 



AM 

Kp 2 ' 3 



H 



Kplin 



7/3 



,7/4 



7/12 



, (KplLl\ 7/3 /Mi 



M 
1/3 /.■ 



M 



r3/4 I Ml 



1/4 



AM 



(105) 



Note that AM /M is the fractional mass annihilated in one free-fall time up to a factor of order unity. 

Consider the fate of a cloud with M/M% < 1, but not <C 1 (we shall estimate how much smaller shortly). Once the 
cloud reaches r > 1, the relativistic products of the annihilation do not escape freely, and deposit a substantial fraction 
of their total energy into kinetic energy of the massive particles. Because the optical depth for the relativistic particles 
is proportional to (E/m) 2 , after a relativistic particle scatters once, thereby losing about half of its original energy 
on average, its optical depth to escape drops by a substantial factor; hence, as long as r is not too large, relativistic 
particles are not trapped. The massive particles that scatter the relativistic particles are heated to semi-relativistic 
kinetic energies i.e. ~ m/2. For these, however, the scattering cross section remains rs ctq] moreover, on each scatter 
off of one of the colder massive particles that dominate the cloud, these hot massive particles lose of order half of their 
kinetic energy. As a result, they continue to slow down, and it is likely that they become trapped inside the cloud. 
The second of equations i|105[l shows that AM/M > Kp 2 ^ 3 for clouds masses <■ Mi at r ^ 1, so the cloud heats up 
to (v 2 ) — T]AM/M, with 77 < 1/2. Thus, we expect the internal energy of the cloud to rise rapidly to i]AM/M as it 
continues to collapse to larger optical depth. 

Since the internal energy of the cloud rises as it annihilates, the annihilation equation is altered to 



dM _ dAM ^ r)pa AM 
dt dt m 



for AM < M we find 



In AM = In AM, + — 
3 



1/2 / \ 1/2' 
P \ (P -, 



vPl/ VPl. 

Taking pi and Mj to correspond to r = 1, we find that AM — » M when the density must increase to 



(106) 



(107) 
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for the cloud to annihilate away. (The last stages of the cloud's collapse could deviate somewhat from free-fall, as the 
pressure builds up to values comparable to pvg, but this phase most likely is only short-lived.) Thus, we estimate 
that the critical mass for a black hole to form is of order 

«„= 2 " M ' - £ fM " m ; (109) 

for example, if m = 50 GeV and gt ~ g* ~ g*s ~ 9*SK ~ 50, we find that M C)P 10(277)M ann or of order 1O 7 (2?7)M0. 
We note that since Mi is of order M ann for s-wave annihilation with (av) ann — cr , the critical mass for black hole 

2/3 2/3 

formation for p-wave annihilation is only smaller by a factor of ~ 7hi[f s /Kp a x m )/2r], not by the factor of ~ KpJ nn 
deduced under the false assumption of adiabatic collapse. 

It is also possible that the CDM particles can annihilate via both s-wave and p-wave channels with p-wave dominant 
in the early Universe but s-wave dominant during the early stages of the cloud collapse. If the s-wave annihilation has 
(ct')ann = e s <7o, then it dominates until AM/M rs e s , whereupon the p-wave channel takes over. In this situation, we 
expect to get a critical mass for black hole formation that is smaller than for s-wavc annihilation with {uv) ann = ctq 
by a logarithmic factor of order (3/2?y) hie" 1 . 



B. Conclusions 



Our calculations show that the collapse of a cloud of annihilating matter to a black hole is only possible for 
masses above a critical mass, M c , that is of order 10 7 — 10 8 M Q for plausible cold dark matter candidates. For 
a radial density profile, the critical mass depends on the parameters of the profile; for example, we found that 
M c /M ann ~ Pa(m for an LL profile with the exponent b |23j . where f3 ann is the density contrast between the center 
and the edge when annihilation becomes important. The conditions for collapse studied here are highly idealized 
since we assumed that the clouds are exactly spherical and cold, at least initially. Internal shearing motions induced 
by tidal perturbations in the early Universe would violate these assumptions. However, it is conceivable that the 
collapse of larger, sheared structures could very rarely lead to the formation of bound fragments similar to what we 
have assumed here mug. In any event, it is likely that we have assumed conditions most favorable for forming 
a black hole from CDM collapse, so our lower bounds on the masses of black holes that can form in this way are 
probably, if anything, underestimates of reality. It is remarkable, although perhaps completely fortuitous, that the 
masses estimated here - which are determined by fundamental constants of Nature - are even close to the sorts of 
black hole masses observed astronomically [l7^. 

When a cloud collapses to a black hole, annihilation releases energy ~ M ann in a time ~ GMbHi where 
Mbh > M ann is the mass of the black hole when collapse first occurs. The peak luminosity resulting from the 
annihilation is therefore ~ M ann /GMBH, and much of the annihilation occurs in a prodigious, relatively short-lived 
burst. Thus the formation of a black hole in the collapse of annihilating CDM would be a spectacular event: for 
s-wave annihilation, approximately M ann <~ 2 x 10 64 [m(GeV)] _1 (a;/-y/g7/15<7* 1 g) ergs would be released in a time of 
order 10 5 [m(GeV))~ 1 (xfy/g^/15g*s){MBH:/M ann ) seconds, with a peak luminosity of 4 x 10 59 (M ann /MBH) ergs per 
second; for m — 50 GeV, the energy released is of order 7 x 10 61 ergs over about 380(MBn/M ann ) seconds. For com- 
parison, note that i(10 11 M Q )(300 km s -1 ) 2 w G(10 11 M Q ) 2 /10 kpc w 10 59 ergs; i.e. the energy released in forming a 
single annihilating CDM black hole is comparable to, if not larger than, the binding energy of an entire galaxy. 

What an observer could detect from such events depends on what annihilation products form, but even if the 
efficiency of generating electromagnetic radiation is < 10~ 8 of the total energy released, they would be as bright 
as extragalactic gamma ray bursts. Any direct production of standard model particles in the annihilation should 
produce a pair fireball with little or no net baryon content, even for extremely low branching ratios for mm — > 
(standard model pairs). For example, if there is no direct production of e ± pairs or quarks, but a fraction f v of the 
annihilation energy goes into neutrinos and antincutrinos of energy to, then we estimate an efficiency of producing e ± 
pairs of about (0.01 — 0.1) f„G 2 F m 2 / (av) ann w (0.001 — 0.01)/^[m(GeV)] 2 (a;/Y / g7/15^s) _1 for s-wave annihilation, and 
perhaps a bit smaller for p-wave annihilation; this would lead to a total energy of order (0.02 — 0.002) x 10 64 / 2 m(GeV) 
ergs 26]. Thus, a tiny branching ratio for mm — » vv would result in the formation of an extremely energetic pair 
fireball, even without any direct pair production via toto — > e + e~ or mm — ► qq: if /„ ~ 10 -5 , we estimate a total pair 
fireball energy of order 10 51 — 10 52 TO.(GeV) ergs. Acceptable branching ratios for direct production of e + e~ or qq pairs 
are far smaller, ~ 10~ 10 . The latest upper bounds on the WIMP-nucleon scalar cross section from CDMS [25| are in 
the range 10" 42 4 - 10 ~ 41 cm 2 w 10~ 15 - lO" 13 - 6 GeV" 2 , compared to (crt;) ann » 8.6 x 10" 9 GeV^ixf^/lbg+s) 
for s-wave annihilation and Co ~ 1-7 x 10~ 7 GeV _2 (x 2 y/g^/150g*s) for p-wave annihilation, so very low branching 
ratios for toto — > (standard model pairs) may be reasonable, if not expected. 
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Astronomical observations could hardly have missed the production of pair fireballs with energy ~ M ann on 
timescales ~ GMbh- Moreover, a blast wave composed of pairs, photons and baryons with a total energy 
that is even a modest fraction of M ann (say 0.001 — 0.01) could deposit enough energy to substantially alter if not 
completely eject the interstellar medium of a galaxy. One possibility is that the branching ratios are < 10~ 10 for 
mm — > e+e~ or qq and < 10~ 5 for mm — » vv, as discussed above. Alternatively, either black holes do not form from 
annihilating dark matter, or else the dark matter is in some other form that does not annihilate, such as axions. 
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